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Abstract 

A simplified model of a redundant power grid is used 
to study integration of fluctuating renewable generation. 
The grid consists of large number of generator and con- 
sumer nodes. The net power consumption is determined 
by the difference between the gross consumption and the 
level of renewable generation. The gross consumption 
is drawn from a narrow distribution representing the 
predictability of aggregated loads, and we consider two 
different distributions representing wind and solar re- 
sources. Each generator is connected to D consumers, 
and redundancy is built in by connecting R < D of these 
consumers to other generators. The lines are switchable 
so that at any instance each consumer is connected to a 
single generator. We explore the capacity of the renew- 
able generation by determining the level of "firm " gen- 
eration capacity that can be displaced for different lev- 
els of redundancy R. We also develop message-passing 
control algorithm for finding switch settings where no 
generator is overloaded. 



1 Introduction 

The design and control of the existing power grid is 
mostly based upon centralized power generation and 
one-way flow of power. The model works quite well 
because inter-hour loads variations are relatively pre- 
dictable allowing reliable centralized scheduling, and 
intra-hour load fluctuations are typically small and can 
be followed efficiently by rapid response regulation re- 
serves. However, the effects of time-variable renewable 
generation are beginning to impact this model, and the 
anticipated growth of wind and photo-voltaic (PV) gen- 
eration at the utility and residential scale could severely 
impact the reliability of electrical power systems (TJ. 

Wind and PV generation are viewed by most util- 
ities and electrical transmission operators as "energy" 



sources, i.e. the electrical system will accept the energy 
that these sources produce when the wind is blowing or 
the sun is shining, but fluctuations in the output of these 
sources are large so that they are not relied upon. This 
is in contrast to a "capacity" source such as fossil, nu- 
clear, or hydroelectric generators whose output can be 
controlled and scheduled to follow the inter-hour load 
variations. When wind and PV are viewed strictly as 
an energy sources, no additional backup generation ca- 
pacity needs to be added. However, the addition of in- 
termittent renewable generation affects the mix and dis- 
patch of other generating and regulating reserves, in- 
creasing the grid's operating cost and the reserve margin 
[3] IU |5] . In essence, the apparent capital and operat- 
ing cost of renewable generation is increased. 

For wind and PV to make up a significant component 
of the generation mix, advances in our understanding, 
planning, and design must be made so that some fraction 
of this new generation can also be counted upon as a 
capacity source |6]|7]|8]. Decreasing the variability, or 
dispersion, of these resources through geographical |9) 
and resource-base diversification or short-term storage 
ifTOlfTT I are possible routes receiving attention. 

We explore an alternative method that utilizes switch- 
able redundant transmission or distribution lines allow- 
ing for reconfiguration of the electrical network to match 
available generation supply to demand without resort- 
ing to load shedding of any type, e.g. line tripping, 
demand response, etc. Specifically, our prime focus is 
on preventing overloads of the firm power generation 
units caused by fluctuations in the demand and inter- 
mittent renewable generation by effectively utilizing an- 
cillary lines. Although the current transmission grid is 
high meshed, we consider a modification where certain 
lines are switchable so that the operational realization of 
the transmission grid is kept radial or tree-like to avoid 
power circulation and certain congestion issues associ- 
ated with meshed grids. Our simplified model is perhaps 
closer to many distribution systems where switchable re- 
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dundant lines of this type already exist and could be ex- 
panded to imitate the type of system we suggest. Note 
also that, as argued in fl2l . switching can play an es- 
sential role in relieving congestion on both transmission 
and distribution levels. 

We explore this method by incrementally adding re- 
newable generation to our simplified electrical model 
and determining the minimum required firm generation 
to ensure that all loads can be satisfied. The reduction 
in the required firm generation as a fraction of the added 
renewable generation provides a measure of its capac- 
ity. We determine the renewable generation capacity for 
different levels of redundancy as well as for different 
amounts of generation variability. 

The simplifications we make to create a tractable 
model are extreme. We completely ignore (a) power 
losses, line impedance, reactive power flow, and tran- 
sient effects llT3ll . essentially modeling electricity deliv- 
ery as an abstract commodity flow fl4l . (b) long dis- 
tance exchange of active power that occurs in a meshed 
transmission-level grid, (c) inhomogeneities and spatio- 
temporal correlations in loads and generation, and (d) 
all economic, pricing and regulation effects fl5l . In 
this very simplified model, we strive to establish the 
fundamental limits or bounds on the capacity of renew- 
able generation and renewable penetration levels beyond 
which no further capacity value is added. 

The material in the paper is organized as follows. In 
Section[2] we introduce our simplified model of the elec- 
trical grid and describe the set of problems addressed in 
this paper. In Section [3] we discuss simple bounds on 
the level of firm generation required for avoiding load 
shedding. Accurate limits of where redundancy allows 
the avoidance of load shedding within our model are 
established in Section [4] We assume that the graph is 
sparse and apply methods of Belief Propagation (BP) to 
create an efficient message-passing algorithm to find an 
acceptable grid configuration among the huge number 
of possibilities. This algorithm, which can also be used 
to control the grid, is discussed in Section [5] Section 
[6] summarizes our approach and proposes future exten- 
sions of the study. 



2 Grid model 

Consider a crude simplification of the power grid with 
M sources/generators each connected to D distinct con- 
sumers, so that the total number of consumers is N = 
MD. Greek/Latin indices will be reserved for genera- 
tors/consumers. For simplicity we assume that each gen- 
erator has a maximum production rate of y, i.e. y = (0 < 
ya < _y|oc = ,M), though inhomogeneities in the 



production can be easily incorporated into our approach. 
The configuration of loads, x — (xt > Q\i — 1, • • • ,N), 
is drawn from an assumed known distribution. As an 
example of the load distribution we will consider the 
joint probability distribution of demand, where each 
load comes from a uniform distribution, centered around 
1 with width A: 



P(x)=Ylp(xi), p(x) 



1/A, |*-l|<A/2 
0, 



otherwise 



(1) 



We generally assume that A <C 1 reflecting the relative 
predictability of electrical loads. 

We also assume that consumers have the capability 
to generate electrical power. We consider two differ- 
ent renewable generation models; one spatially homo- 
geneous and the other inhomogeneous. In the homo- 
geneous model, renewable generation with an installed 
(nameplate) capacity z is present and operational at ev- 
ery node but the output of that generator can fluctuate 
between and z with equal probability, i.e. z = (0 < 
Zi < z\i = ,N) is selected from the distribution 



l/z, 0<z<£ 
0, otherwise. 



(2) 



The large width of 'P(z) compared to V(x) reflects the 
intermittent and highly variable nature of any one re- 
newable generation installation in comparison to the rel- 
atively predictable loads. The homogeneous model is 
representative of wind farms scattered over a relatively 
large area with the distribution ^P(z) describing the vari- 
ability of the wind resource between the different wind 
farms. The form of < P{z) in Eq. [2]is just an example for 
this article. Our methods are amenable to any bounded 
distribution. 

In the inhomogeneous model, we consider that the 
output of the renewable generation at each node is ei- 
ther zero or equal to the nameplate capacity, i.e. 

<B{z)=Y[p{zi), p(zi)=Mz i -z) + (l-fi)8(z i ), 

i 

(3) 

where /, = or 1 . An example of the inhomogeneous 
model is similar to the homogeneous model where we 
replace the wind farms with PV farms and the /, indi- 
cates whether sun at farm i is obscured by clouds or not. 
In our modeling the /; are assigned randomly throughout 
the grid, and we represent the aggregate effect through 
a single variable / = Lfi/N indicating the fraction of 
renewable generation in operation. 

We model the transmission of electrical power from 
generators to consumers as an abstract commodity flow 
ignoring all effects of AC power flow, any correlations 
between fluctuations in load or renewable generation at 
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different consumers, and any economic, pricing and reg- 
ulation effects 11151 . We start with a simple separated 
grid where a group of D consumers is connected to their 
"local" generator, there being M such groups. 

To explore how redundancy affects the capacity value 
of the renewable generation, we consider an intercon- 
nected grid built from the separated one by adding an- 
cillary lines between consumers and generators. For 
our quantitative analysis we choose a random graph 
drawn uniformly at random from the following ensem- 
ble. Of the D consumers connected to each local genera- 
tor, D — R of them are still simply connected to the local 
generator while R of them are also connected to a second 
"remote" generator via an ancillary line, in such a way 
that each firm generator is connected to exactly D + R 
consumers. In total, MR consumers are connected to 
two generators and M(D — R) are connected to one gen- 
erator. See Fig. [T] for an illustration. Although some 
consumers are connected to two generators, the connec- 
tions are assumed to be switchable so that only one is 
utilized so that the resulting "operational" grid is radial 
or tree-like to avoid power circulation and certain con- 
gestion issues associated with meshed grids. 

Without being able to resort to load shedding, we con- 
sider the possibility of redistributing the load via the sys- 
tem of ancillary lines not used under normal (no over- 
load) conditions. Fundamentally, we ask the follow- 
ing question: can an intelligent arrangement of ancillary 
lines among a system of M generators possibly decrease 
the amount of firm generation required system wide by 
giving some consumers the ability to choose a connec- 
tion to a generator other than its local generator. We 
judge the benefit of a certain number of redundant an- 
cillary lines by computing the minimum firm generation 
capacity required to drive the probability of a firm gen- 
erator overload from a finite number to zero in the limit 
of M — > °° and D = 0(1). The notion that a drastic re- 
duction in firm generation is possible stems from basic 
information-theoretic intuition: any finite error proba- 
bility can be reduced to zero via properly introduced re- 
dundancy |16|. 

We describe the connections (active and inactive) be- 
tween loads x and firm generators y via er = (<7, a = 
0,l|{;,a}€ £7), where Q is the bi-partite graph account- 
ing for all generators, consumers and connections. If 
Oia = 1, the connection between load i and generator a 
is active, otherwise it is inactive. We say that a given 
configuration of loads x and renewable generation y is 
satisfiable (SAT) if there exists a matching er such that 
the following set of conditions are simultaneously satis- 
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(a) R = 0,1,2.3. Graph samples. Ancillary connections to 
foreign generators/consumers are shown in color. 

R = l 

A>* Ax Ax A\ 
A>rvfx%^K 




(b) R = 1. Three valid (SAT) configurations (shown in 
black, the rest is in gray) for a sample graph shown in Fig.[TJi. 

Figure 1 : Illustration of the bipartite graph construction 
of our grid with D = 3. 



tied 



(4) 



aedi 



Vae£: £ o ia (xi-Zi) <y, (5) 
ieda 

where di, respectively 3a, stand for all the nodes to 
which i, respectively a, is connected. Eq. (|4]) enforces 
that loads i are only actively connected to a single gen- 
erator, and Eq. |5]l enforces that the sum of all loads, 
net local renewable generation, connected to generator 
a does not exceed the capacity y of that generator. If 
the reverse is true, i.e. there exists no valid er with all 
Eqs. (4]5i satisfied, we say that x is unsatisfiable (UN- 
SAT). 

Let us note here that we can also think about the cur- 
rent loopy power grid as the graph with M(D + R) lines. 
Condition Q is in fact not necessary. However, if there 
are loops in the resulting network of active (switched on) 
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links than the Kichhoff 's circuit laws have to be ensured. 
As long as the resulting network is a tree the power flow 
is equivalent to general commodity flow. Requiring con- 
dition (|4]) can then be regarded as restricting ourselves to 
the configurations of switches for which the Kichhoff s 
circuit laws can be ignored. It is important to notice 
that ignoring condition (HI can only make the satisfiable 
phase larger. We require (|4]) to hold to facilitate our anal- 
ysis. 

First, we aim to solve the decision problem: is the 
given configuration x SAT or UNSAT? Furthermore, if 
x is SAT, we would like to find at least one valid switch- 
ing solution, a. As discussed below, both problems can 
be stated for a given grid or, alternatively, can be consid- 
ered "in average" for ensemble of networks. 

3 Preliminary Considerations 

Let us start discussing some simple bounds on the level 
of firm generation y required to guarantee a SAT config- 
uration. First of all, the total average production has to 
be larger than the total average net consumption. This 
yields for the homogeneous distribution (|2]i of renew- 
ables 

->(1--), (6) 
D ~ 2 

and for the inhomogeneous distribution ([3]) 

| >(!-/£)• (7) 

Another set of bounds can be obtained by considering 
the local structure of the network. If R = 0, we have a 
simple disconnected grid where every generator has D 
consumers. To guarantee that all D consumer loads in 
each group are always SAT without any generator ex- 
ceeding its capacity y and without implementing load 
shedding requires that 

Eq. ([8]) reflects the rare occurrence when every load 
within a local group of D consumers is at its maximum 
1 + A/2 and all renewable generation at the same nodes 
are zero (either because all the zi are zero in the homo- 
geneous model or because all the are zero in the in- 
homogeneous model). Therefore, within this simplified 
model and in the absence of load shedding, the capacity 
value of the renewable generators is zero because they 
do not reduce the required level of firm generation rela- 
tive to the case with no renewable generation, i.e. z = 0. 

This conclusion does not imply that the renewable 
generators have no energy value; they can still be uti- 
lized to offset some of the firm generation y within some 



of the local clusters of a firm generator and D con- 
sumers. However, each local cluster must still have the 
maximum firm-generation capacity y given in Eq. (|8]l to 
avoid all possibility of an overload. 

Allowing one consumer in each cluster an ancillary 
connection to a remote generator (R = 1) yields the same 
result. Consider two local clusters; call them A and B. 
One local load in cluster B is also connected to gener- 
ator A. Additionally, a local load in cluster A is con- 
nected to some other remote generator. We assume that 
each of the D — 1 singly-connected consumers in A and 
B and the one shared consumer all experience the max- 
imum load 1 + A/2 and zero renewable generation. No 
matter how the shared consumer and one other doubly- 
connected consumer in A are switched, either generator 
A or B will always have D connected consumers. There- 
fore, the minimum allowable firm generation capacity y 
is the same as in Eq. ([8j. 

For a level of redundancy of R > 2, the local reason- 
ing is only a bit more complicated. Consider again the 
two clusters A and B. There are order N, 0(N), pairs 
of clusters in the network that share one consumer, and 
0(1) pairs that share two consumers. The probability 
that two clusters share three or more consumer is van- 
ishingly small, o(l), for arbitrarily large but finite R. If 
two clusters share one or two customers, the locally best 
switching that can be implemented will result in at least 
one of the two generators supplying at least D — R+l 
customers. In the rare, but finite probability, event that 
all the singly-connected and shared consumers within 
the A and B clusters are at the maximum load 1 + A/2 
and at the minimum renewable generation of zero, the 
minimum allowable firm generation that still yields a 
SAT state is 

Eq. ([9]) shows that, within the assumptions of our sim- 
plified model, the maximum fraction of firm generation 
that can be displaced by intermittent renewable genera- 
tion is {R—\)/D, and R must be greater than or equal 
to 2 before any capacity benefit is realized. We note that 
this result is achieved solely via switching of ancillary 
connections. Other methods can and should be deployed 
to mitigate fluctuations and allow renewable generation 
to claim a generation capacity factor. 

Eq. (|9]l is obviously not exact as it only sets a bound 
on the maximum amount of firm generation that could 
be displaced by intermittent renewable generation. The 
arguments resulting in Eq. |9} were based on the worst- 
case conditions local to a pair of firm generators, and 
it does not provide any hint as to how much renewable 
generation must be installed to reach this limit. To an- 
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swer this question precisely, we turn to an accurate com- 
putational approach described in the next Section. 



resulting BP equations are 



4 SAT-UNSAT transition 



To solve our model, defined above, we use the cavity 
method and its computational realization via population 
dynamics, introduced in the statistical physics of disor- 
dered systems ifTTlfTSIl and recently adapted to the analy- 
sis of Shannon (phase) transitions in constraint satisfac- 
tion lfT9l and error-correction 12011211 . This method ex- 
plores the famous fact that the Bethe-Peierls |22"1 |231 or 
Belief Propagation ll24l [25 1 (BP) scheme exactly solves 
probabilistic models on graphs without loops (a tree). 
The method allows evaluation of the ensemble averages 
over configurations of allowed discrete switchings on 
the grid that balance load with generation. 

The BP approach to a similar problem, however not 
accounting for a possibility of renewable generation on 
the consumer side, has been explained in [26 1 . Here, we 
provide a brief overview of the technique. In the asymp- 
totic limit of an infinite system, for which N — > °° while 
D and/? are 0(1), the interconnected grid is locally tree- 
like. Therefore, the BP approach for evaluating the gen- 
eralizations of condition (|9]l is expected to be asymptot- 
ically exact (unless the system exhibits a glass transition 
as in |[T8l[T9l . however, we did not observe any signs for 
this breakdown of the asymptotic exactness). We intro- 
duce the following set of marginal probabilities: is 
the probability that generator a is satisfied given that the 
edge (z, a), connecting a with his consumer-neighbor on 
Q, is in the active state, i.e. o !CX = 1. is the prob- 

ability that generator a is satisfied given that the edge 
(z',cc), where (z',0c) € Qi, is inactive, i.e. o ia = 0. %'f^ a 
is the probability that i is satisfied (i.e. it is connected to 
exactly one generator) given that the edge (z',0C), where 
(z, oc) £ Q\, is active, i.e. a, a = 1 . %{p a is the probability 
that z is satisfied (i.e. it is connected to exactly one gen- 
erator) given that the edge (z, a), where (z, a) S Q\, is in- 
active, i.e. Oia = 0. BP relates these marginal probabil- 
ities to each other assuming that the relations are graph 
local, i.e. as if the graph would contain no loops. The 
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(12) 



(13) 



where w; = jc,- + Zi, di \ a is a standard notation for the 
set of generator nodes linked to consumer z, however 
excluding generator a, and similarly 3a \ i stands for 
the set of consumer nodes linked to generator a ex- 
cluding consumer z. The sum over cr^ a \ ia is over all 
values {0,1}I 3 «VI, l n Eqs. l [loffl3| ), Z'^ a and Z a ~ 



are normalizations ensuring that, %' 



1 



Xo 



1 and 



Vo 



1. The 0( ) is the step function enforc- 



ing the generation constraints. It is unity if the argument 
is positive and zero otherwise. The probability for the 
link/edge (z',0C) to be active, stated in terms of the re- 
lated \|/ and %, is 



(14) 



wherel ae3( /7(z,a) = 1. 

The Bethe entropy, defined as the logarithm of the 
number of possible SAT configurations, is 

w = Eiog(z a ) +E lo g( z ') - E lo g( z ' a ) . ( 15 > 

a i (i,a) 

z a = E bMw)I1C < 16 > 

<r 3a ={0,l}l 3 al ieaa ieda 

z'= E^r ! ' n ^ d7) 

ae3i Pe3i'\a 



z — Yl Xl +Y0 Xo 



(18) 



The entropy Seethe is extensive, 0(N), and self- 
averaging, i.e. the distribution of Seethe I s concentrated 
around its mean value with dispersion being o(N) at 
N —>■ oo. 

The fixed point of the BP Eqs. ( T0p"3 i and the cor- 
responding entropy (18 1 can be obtained by solving 
Eqs. ( |T0][T3] > iteratively for a given instance of the prob- 
lem. Repeating these simulations many times for dif- 
ferent instances of the graph and load ensembles, one 
can also calculate the average Bethe entropy. However, 
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the average behavior, corresponding to the limits of the 
infinite graph, can be obtained more efficiently via the 
population dynamics technique. This technique offers 
a computationally efficient sampling from the distribu- 
tion of the marginal probabilities % a on infinite random 
graphs and subsequent evaluation of the average Bethe 
entropy. We will not explain here details of the Popula- 
tion Dynamics approach referring the interested reader 
to DIlEZllSD for further details. 



Implementing the population dynamics method we 
observed three possible outcomes (a) SAT phase: The 
Bethe entropy is positive suggesting that the number of 
SAT configuration (valid redistribution of the demand 
over the generators) is exponential in the system size, 
(b) UNS AT phase, type 1 : The Bethe entropy is nega- 
tive, suggesting that there is almost surely no valid re- 
distribution of demand over the generators, (c) UNSAT 
phase, type 2: A contradiction is encountered in the BP 
equations, formally correspondent to zero values for the 
normalizations in ( fT0|[T3j ). We conclude that the demand 
is incompatible with the graph and respective generator 
assignment. 

We shall also recall here that the SAT-UNSAT transi- 
tion is actually an abrupt transition only in the asymp- 
totic N — * °° limit. In the asymptotic limit N — > °°, the 
SAT phase corresponds to a loss of load probability of 
zero and an expected unserved energy of zero. In the 
UNSAT phase, the loss of load probability is one and 
the expected unserved energy is a finite fraction of the 
total demand. 

Figures [2] [3] and [4] show results of the population dy- 
namics simulations for the homogeneous (wind farm) 
and inhomogeneous (solar) models. Figure [2] depicts 
the required firm generation per customer, y/D, as a 
function of average renewable generation z/2 per cus- 
tomer in the homogeneous model (|2]) for different lev- 
els of redundancy R. The black full line corresponds 
to the condition Q. The horizontal lines to condition 
(|9]l for different values of redundancy R. The lines con- 
necting the symbols indicate the boundaries of the re- 
spective SAT and UNSAT phases obtained via the pop- 
ulation dynamics. These boundaries are established by 
traversing phase space by decreasing y and catching the 
value where the UNSAT conditions are first observed. 
For larger values of z the data points are slightly below 
the boundary line ([9]) due to numerical imprecision (i.e. 
larger population sizes are needed to encounter the rare 
case responsible for condition |9|). 

In the homogeneous model and for the parameter val- 
ues used in Fig. [2] a fraction close to one of the aver- 
age renewable generation can be utilized as capacity for 
level of redundancy R>2. Condition (|9]l combined with 
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Figure 2: Results for the SAT-UNSAT threshold com- 
puted via the population dynamics for the spatially ho- 
mogeneous renewable generation model at different lev- 
els of average renewable generation per consumer node 
. The horizontal lines correspond to condition ([9]). The 
black full line corresponds to the condition (|6j. The col- 
ored lines connecting the data points separate the SAT 
and UNSAT domains lying above and below the lines, 
respectively. Different colors/lines/symbols correspond 
to different values of the redundancy parameter R. Top: 
Example for D — 4 customers per firm generator. Bot- 
tom: D = 3. 

(|6]l gives a rough value of average renewable generation 



R-l 
D 



1- 



R-l 



D 



(19) 



above which adding more renewable generation into the 
grid does not add any extra capacity value. This thresh- 
old increases with the level of redundancy. 

Figure. [3]shows the necessary firm generation per cus- 
tomer for the inhomogeneous model ^ with / = 1/3 
(a third of customers with renewable generation) as a 
function of the amount of the renewable generation ca- 
pability t Again the full black line depicts the global 
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Figure 3: The necessary firm generation per customer 
for the inhomogeneous model with / = 1/3 of cus- 
tomers producing renewables for different levels of re- 
newable generation capability (i.e. nameplate capacity). 
The full black line depicts condition (j7]i and the horizon- 
tal lines condition |9]) for different values of redundancy 
R. The colored lines connecting data points separate the 
SAT and UNSAT domains lying above and below the 
lines, respectively. Different colors/lines/symbols corre- 
spond to different values of redundancy R. 



Figure 4: The necessary firm generation per customer 
for the inhomogeneous model with fixed installed re- 
newable generation capability z = 0.5, 1.0,2.0 as a func- 
tion of the fractional coverage /. The full black horizon- 
tal line this time depicts condition |9]) and the colored 
lines with slopes z represents ([7]). The colored lines, 
connecting the population dynamics data points, sepa- 
rate SAT and UNSAT domains which lay above and be- 
low of the lines, respectively. 



condition |7) and the horizontal lines the local condition 
|9]l for different values of redundancy R. The lines con- 
necting the data points are drawn from the population 
dynamics solution for BP equations. The drop-plateau 
structure of the SAT-UNSAT threshold is related to the 
Independent number of customers which can be served 
by one producer. Qualitatively, the change in behavior 
seen around z ~ 1 is related to customers with renew- 
ables becoming independent of the producers, and an- 
other visible change at z ~ 2 can be explained by the fact 
that consumers with renewables start to supply power re- 
liably to other consumers which do not have renewable 
capability. 

Fig. [4] shows the necessary firm generation per cus- 
tomer computed via the population dynamics technique 
for the inhomogeneous models Q. In this Figure, the 
renewable production capability z is fixed, and the frac- 
tion of consumers / with renewable capability is vari- 
able. The three curves are for different values of t One 
observes that the larger the renewable production z, the 
more one is able to displace firm generation. On the 
other hand, for smaller values of z the fraction of produc- 
tion lost in order to support fluctuations is smaller. The 
drop-plateau structure of the curves can be interpreted. 
The drop at about / = 1/3 for z = 1.0 and z = 2.0 is 
related to the fact that when more than a third of cus- 
tomers has a renewable production then on average one 



customer per producer (D = 3) can become independent 
(for z— 1) or even supply another customer (for z = 2). 
Similar reason applies for the drop at about f = 1/6 for 
z = 2.0, where on average one customer per two produc- 
ers can supply another customer. 

To conclude, this study of the average Bethe entropy 
shows that networks with added ancillary lines are able 
to withstand fluctuations in the renewable generation 
turning this generation into capacity. This effect is am- 
plified with increasing R. The structure of the satisfi- 
able region depends on the distribution of demand and 
renewables production. The drop-plateau structure of 
the curves in the inhomogeneous model suggests that 
even in more realistic systems one might expect a phase 
where adding a few more renewable resources does not 
contribute much additional capacity and a phase where 
adding a few more renewables can, on contrary, add 
more capacity than the amount of additional renewable 
generation, (when the slope of the SAT/UNSAT curve is 
steeper than the slope of lines (|7]l). We did not observe 
this structure in the homogeneous model (Figs. [2]). The 
results of our simplified models are guides to the type 
of generic behavior expected in more realistic models 
which should be investigated in detail to determine ac- 
curate phase diagrams. 
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5 Control Algorithm 

Belief Propagation provides an efficient message- 
passing control algorithm for assigning loads to gener- 
ators. In the spirit of |fl9l , we illustrate that the BP- 
based decimation algorithm is able to find valid config- 
urations within the SAT domain. The decimation algo- 
rithm works as follows: 
DECIMATION: 
1 repeat ([T0p3]i n times 



2 compute the marginals ( 14 1 

3 select the most biased edge and fix its state 

4 cut the fixed edge from the graph 

5 until solution or contradiction is found. 

The BP-based decimation algorithm updates 
Eqs. ( fTOpj} iteratively by passing messages from 
generators to consumers and back. After fixed number 
of steps, the most biased consumer is chosen, the more 
probable value for its consumption is assigned, the 
graph is reduced, and the procedure is repeated. Note 



that updating Eqs. (ToJTTi takes 2 2R steps per message 



thus making the algorithm exponential in R. However, 
building new connections is expensive and R will not be 
very large in realistic cases. The algorithm is quadratic 
in the number of consumers, but linear computational 
time is achievable if a finite fraction of customers is 
decimated at each step, just as in Hl9\ . 




3.3 3.4 3.5 

firm generation, y 



3.7 



Figure 5: Performance of the BP-based decimation al- 
gorithm. The data are for networks with M firm gener- 
ators, n — 5 (five iteration per a cycle of the decimation 
procedure) and R = 2, D = 3, A = 0.67. An average 
over 50 random instances is taken, and the fraction of 
successful runs is plotted against the production cap y. 
The Bethe entropy-based (asymptotic) curve is drawn 
dashed for comparison (the actual value for the curve 
is not related to the success rate), suggesting that in the 
limit of N — > °° valid configurations exist starting from 
y « 3.33. Note that in the separated case (of R = 0), valid 
configurations exist only above 5> > D(l + A/2) = 4. 



The algorithm performance is illustrated in Fig. [5] 
where the percentage of success in the BP-based deci- 
mation is shown. These data average over 50 random 
instances from the R = 2, D = 3 ensemble with a width 
in consumption of A = 0.67, no renewable generation, 
and a variable generation cap. For such a set of parame- 
ters, the separated network would require a firm genera- 
tion capacity of y > 4, whereas with two consumers per 
producer connected to two producers the minimum firm 
generation decreases to y sw 3.33. 

6 Summary and Path Forward 

This manuscript expands upon [26 1 describing a Belief- 
Propagation approach to determining the satisfiability 
of load and generation balance in the absence of load 
shedding for electrical grids including redundant ancil- 
lary connections. Here, we have included intermittent 
renewable generation at consumer nodes and computed 
the minimum level of firm generation required to avoid 
an overload of any the firm generators. The amount of 
firm generation that can be displaced by the combina- 
tion of renewable generation and intelligent switching 
of the ancillary lines equates to the generation capac- 
ity value of the renewable generators. Our results pro- 



vide general guidance about the behavior of renewable 
generation capacity for various levels of redundancy and 
renewable penetration, and we find that the addition of 
ancillary lines greatly increases the capacity value of the 
renewable generators. 

The general approach we pursued in this study is 
based on recent developments in the field of graphical 
models that merge statistical physics, computer science, 
optimization theory and information theory ll20l |2~TI . 
The approach is asymptotically exact on infinite sparse 
graphs and provides an efficient heuristic for graphical 
models on finite sparse graphs. The BP-approach is gen- 
eral and of a dual use: on the one hand BP is a good 
for asymptotic (capacity/phase-transition style) analy- 
sis, and it also generates efficient algorithmic solutions 
for optimization and control of the power grid. 

Note that many generalizations of our analysis are 
very straightforward and can be used directly within the 
framework presented here. This includes implement- 
ing different probability distributions for demands and 
renewable generation (as long as the distribution sup- 
port is bounded) or consideration of non-uniformity for 
both generators and consumers, i.e. varying produc- 
tion caps for either firm or renewable generators and in- 
troducing distinct distributions of demands for different 
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consumers. Also, our simplified network can be easily 
extended in many ways. Our equations are straightfor- 
wardly valid for every locally tree-like random network. 
Another case, which allows very natural and straight- 
forward generalization for all the statements made in 
this manuscript, corresponds to breaking the equiva- 
lence between different edges in the graph and thus as- 
signing nonuniform weights to them. These weights 
may e.g. represent cost of construction, geographical 
length, proxy for losses, cost of exploitation, etc. There 
are many other more realistic extensions of our model 
associated with the description, optimization and con- 
trol of power grids which can benefit from utilizing a 
graphical model approach of the kind discussed in this 
paper, even though actual analysis may prove to be more 
involved. 

Let us emphasize that the BP-based control algorithm 
developed in the manuscript represents a generic ap- 
proach that can be used for efficient switching on any 
power grid, structured or not. This switching algorithm 
is a heuristic which does not rely on any asymptotic as- 
sumption (which is only needed for theoretical analysis 
of the SAT-UNSAT transition), and as such it offers a 
very powerful practical control tool. 
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